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Abstract. We use the calculations derived in a previous paper 
(Mera, Chabrier and Schaeffer, 1997), based on observational 
constraints arising from star counts, microlensing experiments 
and kinematic properties, to determine the amount of dark mat- 
ter under the form of stellar and sub-stellar objects in the dif- 
ferent parts of the Galaxy. This yields the derivation of differ- 
ent mass-models for the Galaxy. In the light of all the afore- 
mentioned constraints, we discuss two models that correspond 
to different conclusions about the nature and the location of the 
Galactic dark matter. In the first model there is a small amount 
of dark matter in the disk, and a large fraction of the dark matter 
in the halo is still undetected and likely to be non-baryonic. The 
second, less conventional model is consistent with entirely, or 
at least predominantly baryonic dark matter, under the form of 
brown dwarfs in the disk and white dwarfs in the dark halo. We 
derive observational predictions for these two models which 
should be verifiable by near future infrared and microlensing 
observations. 

Key words : stars : low-mass, brown dwarfs — The Galaxy 
: stellar content — The Galaxy : halo — Cosmology : dark 
matter 



1. Introduction 

One of the outstanding problems in astrophysics is the nature of 
dark matter in the Universe. Over a wide range of scales the dy- 
namics of gravitationally bound systems reveals the presence of 
a large, dominant amount of dark material. There are in fact two 
dark matter problems. The "large-scale" dark matter that is the 
main contributor to £1 is essentially non-baryonic, if O > 0.1. 
And the "galactic" dark matter, whose presence is inferred from 
the small-scale (velocity dispersion) and large-scale (circular 
rotation velocity) kinematic properties in galaxies. On the other 
hand primordial nucleosynthesis, in particular the concordance 
of the predictions of 4 He, 3 He, 2 H and 7 Li abundances bounds 
f2s, with a lower limit Qb ~ 0.01 h~ 2 (h is the Hubble con- 
stant in unit of 100 km.s -1 Mpc -1 ) which exceeds the den- 
sity of visible baryons (see e.g. Carr, 1994; Copi, Schramm 
and Turner, 1995). Coincidentally, galactic halos span the range 
where dark matter could be entirely baryonic, providing a halo 
radius Rh Ss 50 h^ 1 kpc. 



Within the past few years, microlensing observations in the 
central parts and in the halo of our Galaxy shed new light on the 
galactic dark matter problem by revealing the presence of dark 
objects in the Galaxy and thus of "some" baryonic dark matter 
in the Milky Way. In a previous paper (Mera, Chabrier & Scha- 
effer, 1997; paper I) we have examined the different observa- 
tional constraints arising from star counts at faint magnitudes, 
kinematic properties and microlensing experiments to derive 
mass functions in the Galactic disk and halo and to determine 
as accurately as possible the stellar and substellar contributions 
to the Galactic mass budget. 

In this paper, we examine different mass-models for the 
Galaxy compatible with these determinations and thus with all 
presently available observational constraints. In §2, we first dis- 
cuss the nature of the possible stellar content of the halo. In 
§3, we examine what we call the conventional scenario, with 
a standard massive halo and a substantial fraction of unidenti- 
fied, most likely non-baryonic matter. In §4, we devote some 
attention to the maximal disk scenario, whose interest has been 
renewed in the last years. Section 5 is devoted to the derivation 
of a less conventional model, which determines the limits of an 
entirely-baryonic dark matter model for the Galaxy. Summary 
and conclusions are outlined in §6. We propose observational 
predictions characteristic of these models which should be in 
reach of near future infrared and microlensing projects. 

2. Which objects in the dark halo ? 

The first analysis of the Hubble Space Telescope (HST) deep 
star counts showed that M-dwarfs provide at most ~ 10 — 20% 
of the halo dynamical mass (Bahcall et al., 1994; Santiago, 
Gilmore & Elson, 1996). A more detailed analysis, consis- 
tent with the observed number of high-velocity subdwarfs in 
the solar neighborhood (Dahn et al, 1995) places much more 
stringent limits with a main sequence star contribution to the 
halo mass smaller than 1% (Mera, Chabrier & Schaeffer, 1996; 
Graff & Freese, 1996; Chabrier & Mera, 1997a; Mera et al., 
1997, paper I). Therefore M-dwarfs are definitely excluded as a 
significant halo population and can be responsible for less than 
0.1 of the 6 to 8 events]] observed towards the LMC (see paper 

1 It has been suggested by Gould et al. (1997) that one of these 
events might be due to a disk M-dwarf. 
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I). This raises a severe problem for the observed microlensing 
events towards the LMC. The nature of the objects responsible 
for these events is still under debate, and represents a central 
issue for any model of the Galaxy. 

The optical depth inferred from the LMC observations is 
r Q bs = 2.2 ± 1 x 1CP 7 , and the average mass for the lens is 
(y/m) 2 s» 0.5 M (Alcock et al., 1996b). This mass is deter- 
mined with a standard, spherically symmetric, halo model (see 
paper I §4.3). The velocity dispersion decreases with oblate- 
ness, leading to a lower average mass. Within the range of 
flattened Evans models (Evans, 1993; 1994), the characteris- 
tic mass of lenses remains larger than 0. 1 M Q (see Paper I 
§4.3). Only with a very flattened halo (axis ratio of the order 
of 1/10), the velocity dispersion is low enough for the aver- 
age mass to lie under the hydrogen burning limit. Moreover, 
such a disk-like halo must have a relatively small asymmetric 
drift, which yields an even lower ( 1 / v± ) ~ 1 and an average lens 
mass (\/m) 2 as 0.08M Q . But the spheroid MF normalization 
at 0.1 M Q (Eq. (16) of paper I) is reasonably well constrained. 
Chabrier and Mera (1997a) have shown that the local spheroid 
subdwarfs are not related to the dark matter halo, yielding an 
even more severe constraint on the dark halo MF, whose max- 
imum normalization is ~ 1% of the spheroid one at 0.1 M Q . 
Moreover, the extrapolation of the spheroid MF into the brown 
dwarf domain down to a minimum mass corresponding to an 
average mass of 0.08 M Q yields a mass density of less than 
2% of the dark halo density. Clearly, the substellar part of that 
MF would have to be a Dirac delta-function at 0.08 M Q for the 
whole observed LMC optical depth to be due to brown dwarfs 
and these brown dwarfs should be visible in the Hubble Deep 
Field. In any case a highy flattened halo raises severe problems. 
Features like the Magellanic Stream, for example, can not be 
accounted for in such a model. Reasonable values for the halo 
oblateness are more in the range q ~ 0.6 (e.g. Sackett et al., 
1994), so that brown dwarfs seem to be definitely excluded as 
halo dark matter candidates, unless they are strongly inhomo- 
geneously distributed in dark clusters (Kerins, 1997). 

The expected number of microlensing events in the LMC 
itself is: 



Nlmc 



25 km.s" 1 V 250 VC 36OM pc- 2 



(1) 



where a, E and h denote the velocity dispersion, the surface 
density and the depth of the LMC []. Reasonable values for 
these parameters (Gould, 1995) yield Nlmc ~ L which is 
consistent with the analysis of the binary event in the recent 
MACHO results (Bennett et al., 1997). 

The only remaining (though hardly satisfying) solution to 
explain the observed events toward the LMC is massive stellar 
remnants, white dwarfs (WD) or neutron stars (NS), as first 
suggested by the MACHO collaboration. Halo NS could be 
disk NS whose high-proper motion exceeds the disk escape 
velocity. However, the total number of NS in the Galatic disk 

2 Note that these variables are not independent from each other. 



amounts to ~ 10 9 , with a NS mass strongly peaked around 
~ 1.4M Q . This yields a negligible (<< 1%) contribution to 
the observed optical depth. Moreover, NS are the end product 
of type II supernovae, whose expected number in the Galac- 
tic halo is severely constrained by the observed metallicity. 
In any case NS's and stellar black holes are likely to con- 
tribute considerably less than WDs to the dark mass fraction. 
If the ratio of WDs to low-mass stars in the halo is assumed 
to be the same as the one observed in the disk, i.e. ~ 1/100 
(NwD dtsk ~ 3 x lO^pc" 3 , Liebert, Dahn & Monet, 1988) 
and if we assume the same amount of mass in NS - a very up- 
per limit -, we get a maximum contribution of baryonic stellar 
remnants to the Galactic dynamic mass of ~ 5%. This cor- 
responds to less than 1 event towards the LMC. Clearly, for 
WDs to be responsible for the dark events towards the LMC, 
one must advocate very particular conditions in the primordial 
halo. 

Recently, Chabrier, Segretain & Mera (1996) and Adams 
& Laughlin (1996) have examined in detail the possibility for 
WDs to provide a significant halo population. These authors 
have shown that it is indeed possible for these objects to be 
responsible for the observed MACHO events, thus providing 
~ 40% of the missing mass, under two necessary conditions, 
namely i) a halo age substantially older than the disk and a very 
peculiar halo initial mass function (IMF) strongly peaked (or 
cut) around ~ 1.5 — 2M , completely different from the one 
determined in globular clusters and in the spheroid (Chabrier 
& Mera, 1997a) or in the Galactic disk (Mera, Chabrier & 
Baraffe, 1996; Mera et al., 1997, paper I). This scenario also 
implies that a substantial fraction of the helium produced in the 
WD progenitors has been blown away from the Galaxy along 
its evolution. This is supported by the recent suggestion of the 
presence of hot, metal-rich gas in the Local Group (Suto et al., 
1996; Fields et al., 1997), as observed in other groups of galax- 
ies, although this scenario has been questioned recently (Gib- 
son & Mould, 1997). Although the WD hypothesis seems very 
speculative and raises important questions about the nucleosyn- 
thesis of primordial intermediate-mass stars (Gibson & Mould, 
1997) and requires unambigous confirmation of microlensing 
experiments in the Galactic halo, it cannot be excluded a pri- 
ori. As shown by Chabrier et al. (1996), it is consistent with 
constraints arising from the expected radiation signature of the 
progenitors at large redshift (Chariot & Silk, 1995) and from 
the observed (and unobserved !) density of high velocity white 
dwarfs in the solar neighborhood (Liebert et al., 1988). Cal- 
culations of the expected discovery functions and star counts 
at faint magnitudes for various halo white dwarf luminosity 
functions are under way (Chabrier and Mera, 1997b), that will 
be confronted to near-future high-proper-motion surveys. This 
will soon confirm or rule out the hypothesis of a significant halo 
WD population. 

Interestingly enough, the WD scenario also bears important 
consequences on our understanding of the formation of Type la 
supernovae. If SNI form from accreting WDs, a high-density 
WD halo would imply 100 times more supernovae in the halo, 
with M WD w 0.40 x 10 12 M Q within 50 kpc, than in the disk, 
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where M WD w 0.2 x 10 10 M Q . Such a high rate of SNI in the 
halo of our Galaxy and neighbor galaxies is clearly excluded 
by the observations. In case WD's do indeed represent a sig- 
nificant halo population, this may imply that SNIa do not form 
from accreting WDs, as suggested by several authors (Benz et 
al. 1990; Mochkovitch & Livio, 1989, 1990; Nomoto & Iben, 
1985; Saio & Nomoto, 1985). The final outcome of merging 
WDs might rather be a single, massive white dwarf (Benz et al. 
1990; Segretain, Chabrier & Mochkovitch, 1997). Another pos- 
sibility, in case direct observations confirm the WD hypothesis, 
could be that their mass function is sufficiently narrow to pre- 
vent any instability due to accretion, or that they are old enough 
so that mass exchanges in binary systems no longer occur. 

Although other possibilities may be considered for the mi- 
crolensing objects, such as, for instance, new (unknown) vari- 
able stars, hypothetical novae eruptions (Delia Valle and Livio, 
1996) or primordial black holes, white dwarfs remain the least 
unlikely solution. The models for the Galaxy can thus be split 
into two types: the first one postulates that the dark halo is made 
of matter of unknown nature and does not necessarily require 
the halo WD scenario (§3), the second type considers the pos- 
sibility for WD-like objects to explain all, or at least most of 
the halo dark matter (§5). 



3. Standard massive halo model 

The determination of the different contributions to the dark 
matter density in the solar neighborhood can be understood by 
the following simplified analytical approach. 

For a rotation velocity independent of the galactocentric 
distance, the Poisson equation implies a density distribution for 
the dark halo, assuming it is the only contributor to the dynam- 
ics : 



Ph(r) 



1 



4ttG' 



(2) 



where t> rot ~ 220 km.s 1 is the asymptotic rotation velocity 
(see paper I). The maximal halo solar density corresponds to: 



(max) 
PhQ 



1 



AttG R% 



0.012 M Q pc~ 



(3) 



where R@ = 8.5 kpc is the galactocentric position of the 
Sun. The disk, however, contributes to the rotation velocity in 
the solar neighbourhood v®. This implies a slightly reduced 
local halo density as compared to the maximal halo ([}]). At 
large distances, the disk exponential density-profile yields a 
much smaller correction. Since there is no hint for a signifi- 
cant drop in the observed rotation curve up to at least 70 kpc, 
this compensation is usually taken into account by introducing 
a screened halo density-profile: 



Ph(r) 



1 



attg m 



(4) 



which reduces the halo density in the solar neighborhood to 
(for R c = 5 kpc): 



(max) _ V rot 

PhQ ~ AirG EZ 



1 



R 2 



= 9.2 x lO" 3 M pc" 



(5) 



The rotation velocity v Q is thus related to the halo local vol- 
ume density pn & and to the disk oc exp(— r/Rd) (where Rd is 
the disk scale length) local surface density E , approximately 
by (see e.g. Binney & Tremaine, 1987): 



2ttG 
where 



2 /(-R c / 'Ro)R%ph c . 



M, 



bule 



-disk 



2ttG 2-kRc 



(6) 



o 



f(x) = [l-a;Arctan(l/x)][l + x 2 ] 

v 2 disk /2nG = 2R d y 2 [I (y)K (y) - h (y)K x (y)} S e^/^ 

y = R /2R d , 

I n and K n are Bessel functions, and Afbuigc is the mass of the 
central bulge. 

On the other hand, the contributions of the disk and the halo 
to the vertical acceleration are given schematically by (assum- 
ing z to be small enough for A 2 — B 2 as 0, where A and B are 
the Oort constants): 



1 



K{z) 

2ttG 2ttG dz 



(7) 



The standard halo model corresponds to the self-consistent 
solution of equations (Q) and (^ which determines R c and E 
from the measured u ro t and -jq - For v ro t — 220 km.s -1 , 
^§ = 71M Q pc- 2 at z = 1.1 kpc, R Q = 8.5 kpc and 
R d = 3.5 kpc, we get R c = 4.1 kpc and S = 49 M Q pc -2 , as 
Kuijken and Gilmore (1991). The exact determination of the re- 
spective contributions to the gravitational acceleration (0) due 
to the disk and halo is in fact more complicated and requires de- 
tailed (model-dependent) numerical computations (see e.g. Bi- 
enayme, Robin & Creze, 1987), but the present approach yields 
reasonnably accurate values for the dynamical constraints to 
be fulfilled by different disk/halo models. It is clear that all 
the values of the disk surface density derived from the mea- 
surement of the vertical acceleration in the literature were done 
within the standard halo model, and thus include implicitly this 
assumption. Any Galactic mass-model using a different (non 
maximal) halo will thus imply values for the observed surface 
density larger than the ones currently given in the litterature 
(since the halo contribution will be smaller than the present 
standard one). This important issue will be examined in §4 and 
5. 

We argued in Paper I that the most recent determination 
of the surface density by Flynn and Fuchs (1994), E Q = 
52 ± 13M pc~ 2 , which we use in the present calculations, 
is consistent at the la level with all previous determinations 
and is nearly identical to the uncertainty-weighted average 
S = 51 ± 6M pc -2 of these different but consistent mea- 
surements. The density under the form of observed objects, 
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determined in paper I (Eq. 15), £q = 43 ± 5M Q pc~ 2 , is 
compatible at the ler level with the afore-mentioned dynami- 
cal determination, leaving limited possibility for a brown dwarf 
component Sbd = 9 ± 14 M Q pc~ 2 . If, however, the ler up- 
per limit is retained, then the disk may contain a brown dwarf 
density as high as £bd q — 22Mqpc~ 2 , comparable to the 
main sequence star contribution, still compatible with the mi- 
crolensing observations towards the bulge, as shown in section 
4. 1 .2 of paper I. 

The final parameters for the mass distribution in the Galaxy 
in the conventional model are thus: 

• an exponential disk with a scale length Rd = 3.5 kpc, 
the surface density and brown dwarf contributions mentioned 
above, and a mass M d = 27r£ Qj R 2 , exp(i? Q /i? d ) w 4.6 x 10 10 
M , 

• a central bar (bulge), whose maximum mass is deter- 
mined below, which contributes to the microlensing optical 
depth with Eq. (25) of Paper I, 

• a spheroid with a (flattened) De Vaucouleurs profile and 
a mass density deduced from the mass function given by Eq. 
(16) of paper I. This yields a mass M sph = 1.3 x 1O 1O M , 
which contributes ~ 1% of the local dark matter density and 
—0.1% of the disk total density (Chabrier & Mera, 1997a), and 
an optical depth r sp h — 5 x 10 -9 towards the LMC. 

• a screened standard halo described by Eq. (Q), which cor- 
responds to a halo mass (see §2 of Paper I) Mh = v^ ot r/G w 
10 12 (r/100kpc) M Q . As shown in paper I (see also Chabrier 
& Mera, 1997a) the contribution of main sequence stars and 
brown dwarfs to the halo mass budget is negligible, if not zero. 

In this model, one has to invoke a substantial amount of 
dark matter of unknown origin. 

This model gives the rotation curve shown in figure § We 
see immediately that the contribution of the central bar to the 
velocity at 1 kpc from the Galactic centre cannot exceed the 
observed value of 230 km.s - . This sets the mass Afbuigc ~ 
1.2 x 10 10 Mq within 1 kpc, under the assumption of spherical 
symmetry. The presence of a central bar changes this result to 
M bu ige - 2 x 10 10 M Q (Zhao et al., 1995). 

In summary, the standard halo model is consistent with all 
observed properties of the Galaxy, at the price of postulating 
the existence of a sizeable fraction of non-baryonic matter in 
the halo (that is needed anyway at larger scales if the density 
parameter of the Universe fl is larger than 0.1). To explain the 
observed microlensing events towards the LMC, the only (un- 
comfortable) possibility is to invoke a substantial amount of 
white dwarfs in the halo, which contribute to — 40% of the dy- 
namical mass, providing very particular hypothesis on the age 
and the initial mass function of the halo, as discussed in the 
previous section. Since, as mentioned above, brown dwarfs and 
low-mass stars seem to be excluded as a significant halo popu- 
lation, this model implies that — 60% of the Galactic missing 
mass is in the form of undetected non-baryonic matter. This 
fraction may be as high as 90% if the lower limit for the bary- 
onic fraction in the halo inferred from the number of microlens- 
ing events is retained. Note, however, that this model does not 
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Fig. 1. Rotation curve of the standard model (see text). The 
open circles correspond to the averaged observed values (from 
Vallee, 1994), the dashed line is the disk contribution, the dash- 
doted line the halo contribution, and the solid line is the total 
model rotation curve. 



need a white dwarf component in the halo. Would the future 
number of microlensing events towards the LMC decrease, or 
even vanish, this model remains valid. 



4. Maximal disk 

It has been proposed recently on various grounds (e.g. Pfen- 
niger, Combes & Martinet, 1994) that, contrary to the original 
suggestion of Ostriker and Peebles (1973), the Galaxy has no 
massive dark halo. In this case, all dark matter is distributed in 
a strongly flattened, disk-like structure. 

It is well known that an exponential disk gives rise to a 
decreasing rotation curve, incompatible with the observed flat 
rotation curve (see e.g. Binney & Tremaine, 1987). The inverse 
problem, which consists of deriving a density profile from an 
observed rotation curve, can be solved analytically for the pe- 
culiar case of a constant rotation velocity (Mestel, 1963). For 
an infinite disk, the surface density is proportional to the radial 
distance 1 /R: 



R 



with 



208 i 



Wrot 



220 km.s" 



M pcf 



(8) 



(9) 



This value is much larger than the one inferred from the 
observed vertical acceleration (see Paper I). 

As mentioned above, Eq. (||) holds for an infinite disk. The 
finite-size yields a smaller value for S Q for the same rotation 
curve. In this case, the density profile is very close to an expo- 
nential, S(i?) Y, e- R / Rd , with a scale length R d « 1/3 R G , 
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where Rq is the Galactic (disk) radius (Mera, Mizony & Bail- 
Ion, 1997). Only the very central and external parts of such a 
maximal disk depart from the exponential profile, as expected. 
In the centre, this departure corresponds to the bulge whereas 
at the edge, the density must decrease to zero. Thus, for a trun- 
cated disk, a flat rotation curve does not correspond to a 1/R 
surface density but rather to an exponential disk profile plus a 
central bulge (see Mera, Mizony & Baillon, 1997). The scale 
length of this exponential is Rd ~ 10 kpc for spiral galaxies 
similar to the Milky Way, which corresponds to an exponen- 
tial variation of the mass-to-light ratio M/L in the Galactic 
disk, and an increasing dark matter/visible matter ratio from 
the centre to the edge of the disk. As demonstrated by Mera 
et al. (1997b), if all the mass is distributed in a disk, the local 
surface density is fixed essentially by the disk radius Rq- For 
the rotation curve to be flat up to at least R& = 30 kpc, we get: 



£ > 19OM pc" 



(10) 



As stressed in §2, the "observed" dynamical densities are 
derived by measuring K(z) at some altitude z, corrected from 
the standard halo contribution to obtain the disk density. If the 
halo local density is reduced, the disk surface density implied 
by a given vertical acceleration is different (Eq. For the 
standard halo (§3), the halo correction yields : 



K({z)) 
2ttG 



20 M q pc" 



1 kpc 



(11) 



-2 (z) 



where (z) is the average height of the sample. 

In case there is no halo (maximal disk), 20 M pc~ M . |)( 
must be added to the £ values quoted by the various au- 
thors. With our estimate of (z) = 1 kpc, 1.1 kpc, 0.5 kpc and 
0.5 kpc for the Bienayme et al (1987), Kuijken and Guilmore 
(1991), Bahcall et al (1992) and Flynn and Fuchs (1994) sur- 
veys, respectively, the observed values become S = 70 ± 10, 
70±9, 95±25, 62±13M pc~ 2 , respectively, with an average 
70 ± 6M pc -2 []. Therefore, the surface density ( |ic| ) needed 
to comply with the rotation curve is many standard deviations 
off the observed value. 

One may argue that the vertical distribution of the gas in- 
creases with Galactocentric distances, as expected if the gas 
is supported by rotation (Pfenniger et al, 1994). This would 
yield a dark "corona" beyond the edge of the observed opti- 
cal disk. However, this is not a satisfactory solution, since, as 
can be shown easily, adding a dark corona of matter outside the 
disk will lower the rotation velocity inside the Galactocentric 
radius of this corona, so that the disk surface density must be 
increased accordingly in order to reproduce the observed rota- 
tion velocity, to a value close to ([}]). 

In sum, the total mass of the disk in this model is Md « 
10 11 x i?o/10 kpc. It seems unrealistic that the disk extends 



3 As mentioned in Paper I, the Bahcall et al. (1992) and Flynn & 
Fuchs (1994) observations are not completely independent but the 
weight of the Bahcall et al. value is small and droping it does not 
affect the result. 



beyond 40 kpc for stability reasons so that such a description 
of our Galaxy would not be consistent with the velocity fields 
observed up to 60 kpc (see §2 of Paper I). It similarly implies 
that the mass of the Milky Way is significantly smaller than the 
required 2.4 ± 0.8 x 10 12 M deduced from the dynamics of 
the Local Group (Peebles, 1994; see §2 of Paper I). Or it would 
impose that almost all the mass of the Local Group (~ 4x 10 12 
M ) is in Andromeda. The dynamics of our Galaxy would be 
explained without the presence of non-baryonic dark matter at 
the price of an unacceptably large amount of such material in 
Andromeda. 

The only way out is to require the dark matter to be dis- 
tributed over a scale height z > 3 kpc. In this case, the maxi- 
mal disk is similar to a highly flattened halo, a solution appar- 
ently not supported by observations, as mentioned in §2. If the 
sought baryonic matter is in the form of dark, molecular hydro- 
gen clouds (Pfenniger et al, 1994; Gerhard & Silk, 1996), this 
raises a severe problem, for the observed distribution of atomic 
hydrogen (Dickey and Lockman, 1990) is well confined within 
200 pc from the Galactic plane. This implies that molecular and 
atomic hydrogen gas must have a completely different (a factor 
10) Galactic scale height. 

In summary, the maximal disk or highly flattened halo mod- 
els seem to be excluded by kinematic observations and a 3- 
dimensional dark halo contribution is necessary, which is re- 
sponsible for the constant (non-keplerian) rotation curve be- 
yond the disk truncation radius. 



5. Modified halo model 

As shown in Eq. (Q), the halo distribution is determined by two 
parameters: its size, fixed by the total Galactic mass, and its 
core radius R c . This latter is not physically motivated but is de- 
termined by the condition that the visible plus dark matter dis- 
tributions give rise to the measured rotation curve, which yields 
the local halo normalization. The standard halo model (§3) as- 
sumes by definition that the dark matter is mostly in the dark 
halo, so that the local halo density is ph e = 9 x 10~ 3 M pc~ 3 
and R c = 5 ± 1 kpc (Eq. g). 

We modify the standard model by putting as much baryonic 
mass in the disk as allowed by the microlensing observations 
and the star counts, from the analysis conducted in Paper I. 
The aim is to derive the maximum limit of a predominantly 
baryonic model for the Galaxy compatible with all presently 
available observations. 

We first note that there is observational evidence that the 
optical disk is truncated at a galactocentric distance R w 14 
kpc (Robin, Creze & Mohan, 1992; Ruphy et al., 1996). As 
mentioned in section §4, such a truncated disk plus a central 
bulge yield a flat rotation curve up to the limit of the disk, with 
a scale length Rd ~ 3.3 kpc, in agreement with observations 
(1.8 < R d < 5 kpc, see Kent, Dame & Fazio, 1991) and S w 
12OM pc~ 2 (Mera et al., 1997b). Since a maximal disk is 
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excluded (§4), we must include a halo contribution given by 
the following flattened distribution (in cylindrical coordinates): 



ph{R,z) 
with 

Ph @ = f X 



R r . 



R 2 + z 2 /q 2 + R 2 



AttG R 2 



Rq q Arccosq ' 



(12) 



(13) 



where / is the normalization of the halo, and q ~ 0.6. Straight- 
forward calculations yield the following expression for the con- 
tribution of such a halo to the rotation curve: 



2 

rot 



Rc6 Arctan^e 2 + R 2 /R 2 C 



R 2 



Arcsin e 



(14) 



q 2 . For R >> R c , the asymptotic rotation 



where e = y/l 

velocity is \ffv ro t, with w rot « 220 krn.s -1 . 

The core radius R c plays in fact the same role as the 
normalization factor /, as shown by eqn. (|l2|). Halo mod- 
els (/ = 0.35; R c = 5 kpc) and (/ = 1;R C = 15 kpc) 
yield the same normalization in the solar neighborhood ph Q = 
3 x 10~ 3 M pc~ 3 , i.e. ~ 30% of the standard model local 
density, but the second one yields the correct rotation velocity 
at large distances (R <; 20 kpc). A truncated disk plus a dark 
halo with a large core-radius, R c ~ 15 kpc is thus a solution 
worth considering. This model is consistent with the microlens- 
ing experiments towards the LMC which yield a dark matter 
contribution ~ 40% of the standard halo: the microlensing op- 
tical depth toward the LMC derived from the afore-mentioned 
density is r = 3 x 10~ 7 , in agreement with the observed value. 

We have substracted the contribution of this modified halo 
from the observed rotation curve, taken from Vallee (1994), up 
to ~ 14 kpc. Figure 2 shows the inferred density profile for 
the disk. The agreement between theory and observation yields 
a scale length Rd — 3.2 kpc and a disk mean surface density 
E Q = 87M Q pc- 2 . 

Instead of using the numerical density profile derived from 
the inversion technique of Mera et al. (1997b), we will use the 
more simple exponential density law. Then the observed rota- 
tion curve is well reproduced with the following disk density: 



E d (R) = (75 Mq pc- 2 ) x 



e R e /R d e -R/R d 



Figure 3 shows the total (disk + halo) rotation curve up to 
60 kpc. The disk dominates the rotation curve up to its edge, 
i.e. ~ 14 kpc, whereas the dark halo becomes dominant beyond 
this limit and reaches the correct asymptotic value 220 km.s -1 . 
The rotation curve of the model departs significantly from the 
observations for R > 12 kpc, but the uncertainty begins to 
be large at these radii (about 20%). In particular, the last two 
points are derived from 5 objects only, each with large errors. 
Therefore, we can consider the model as consistent with the 
observed rotation curve. 

The most interesting property of this model is that it is con- 
sistent with a fully baryonic dark matter solution, most of the 




4 6 8 10 

Distance from the center in kpc 



14 



Fig. 2. Surface density profile of the Galactic disk inferred from 
the observed rotation curve (Vallee, 1994), from which the 
modified-halo (R c = 15 kpc) contribution has been removed 
(see text). The dashed line is an exponential profile of scale 
length 3.2 kpc, and solar density 87 M pc -2 . 
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Fig. 3. Rotation curve of the modified halo model. The disk 
contribution corresponds to the dashed line, while the halo con- 
tribution corresponds to the dash-dotted line. The solid line 
is the total rotation curve, with the observed rotation curve 
(Vallee, 1994) for R < 14 kpc. Note that the observational 
uncertainties become very large beyond 10 kpc (see Fich & 
Tremaine, 1991). 



matter being in the dark halo, and some of it in the disk. The 
disk dark population consists of brown dwarfs (see §4.1.2 of 
Paper I), whereas for the halo, the most likely candidates are 
white dwarfs, as discussed in §2. As discussed at length in §3.1 
and §4.1 of Paper I, a substantial population of disk brown 
dwarfs is consistent with present-day microlensing observa- 
tions toward the Galactic center. It is also consistent with the 
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extrapolation into the sub-stellar domain of the M-dwarf mass 
function obtained either from the nearby luminosity function 
(Kroupa, 1995) or from the very faint end of the HST lumi- 
nosity function (Gould et al., 1997), if the apparent upturn in 
this region is real. Such mass functions dN/dm oc m~ a have 
a slope a <; 2 and a normalization near the bottom of the main 
sequence dN/dm(Q.l M©) ~ IMq^c -3 (see §3.1 of Paper 
I). 

A rising MF near the brown dwarf domain is supported by 
the recent study of the distribution of secondary masses of 20 F 
and G dwarf binaries by Mazeh, Latham & Stefanik (1996), al- 
though poor statistics (3 detections) preclude the derivation of 
a precise MF. On the other hand, a shallower MF or a substan- 
tially smaller normalization near the bottom of the MS will rule 
out the possibility of a substantial brown dwarf mass-fraction in 
the Galactic disk. Therefore nailing down the issue on the very 
end of the disk MF bears essential consequences for Galactic 
modelling. This should be done shortly with the DENIS whole- 
sky survey. 

As shown above, the main problem of such a modified-halo 
model is the large inferred surface density £© ~ 75 M Q pc~ 2 . 
For this model, the halo correction is only 



K({z)) 
2ttG 



6Mq pc~ 



(z) 



1 kpc 



(15) 



-2 (z) 



the 



so we have to increase by 14 (20 minus 6) M Q pc - ^ 

observed dynamical values. This yields £© = 63 ± 10, 62 ± 
9, 91 ± 25, 58 ± 13 M p C - 2 for Bienayme et al (1987), 
Kuijken and Gilmore (1989), Bahcall et al (1992), Flynn and 
Fuchs (1994) determinations, respectively, with the average 
S Q = 63 ± 6M Q pc -2 . Since the detected surface density 
is E 0t , ia = 43 ± 5M pc" 2 (Eq. 15 of Paper I), the brown 
dwarf surface density in the solar neighborhood is then Sbd = 
15 ± 14, still compatible with microlensing observations (see 
§4.1.2 of Paper I). 

The axisymmetric mean density at the solar radius E© w 
75 M Q pc~ 2 obtained in this model is larger by ~ 1.3 standard 
deviations than the afore-mentioned Flynn & Fuchs (1994) re- 
cent determination in the solar neighborhood. The difference, 
however, may be attributed to the fact that the Sun is located 
in-between two spiral arms, namely ~ 1.5 kpc from the Sagit- 
tarius arm and ~ 2 kpc from the Perseus arm (King, 1989). 
Therefore the solar surface density might be slightly smaller 
than the averaged surface density along the circle of radius i?©. 
This is supported by the observations of external galaxies and 
by numerical simulations, which give about a factor 2 for the 
maximal/minimal density in the disk, expected to translate into 
a smaller factor for the vertical acceleration K{z) (see e.g. Rix 
and Zaritsky, 1995). The value at 1.1 kpc of K z /2nG derived 
by Kuijken and Gilmore (1991) is 71 ± 6 M Q pc~ 2 , whereas 
for the present heavy disk model K z /2nG — 81 M© pc~ 2 , 
corresponding to a ratio of 1.14, which could be accounted for 
by the spiral structure. 

The average density we need is ~ 1.3 times larger than the 
mean solar value, and only 5% times larger than the 1-a up- 



per limit of the Flynn and Fuchs (1994) value. Although the 
relevance of such values deserves a detailed study, they do not 
seem to be unrealistic. 



6. Summary and conclusions 

In the present paper, we have derived consistent mass-models 
for the Galaxy, disk(s), bulge, spheroid and dark halo. The 
present calculations differ substantially from previous ap- 
proaches which relied on microlensing analysis only (e.g. Al- 
cock et al., 1996b) or which were based on the minimization 
of models with a large number of free parameters (e.g. Gates, 
Gyuk and Turner, 1996). We proceed differently by first de- 
termining as accurately as possible the main parameters enter- 
ing the Galactic modelling from all observational constraints 
presently available, i.e. star counts and luminosity functions of 
metal-rich and metal-poor M-dwarfs at faint magnitudes, in- 
cluding the most recent Hubble Deep Field data, microlens- 
ing observations towards the LMC and the central part of the 
Galaxy, including not only the optical depths but the much 
more constraining time-distributions of the events, and large- 
scale (circular rotation velocity) and small-scale (local gravi- 
tational acceleration) kinematic properties. The microlensing 
analysis is derived not only in term of optical depth but in- 
cludes the complete time-distributions of the events. These cal- 
culations include previous results wich yield the determina- 
tion of the mass-functions in the Galactic disk (Mera, Chabrier 
& Baraffe, 1996) and spheroid (Mera, Chabrier & Schaeffer, 
1996; Chabrier & Mera, 1997a) and of the contribution of stel- 
lar and substellar objects to the Galactic disk and halo mass 
budgets (Mera, Chabrier & Schaeffer, 1997; Paper I; Chabrier 
& Mera, 1997a). This yields a consistent determination of the 
different parameters. 

We show that a maximal disk model is clearly excluded 
by the observations, in particular by the observed local surface 
density. An extension of this model, where dark matter is in the 
form of giant molecular clouds located at the very periphery of 
the Galactic disk (Pfenniger et al, 1994) raises severe problems 
: i) it implies that molecular hydrogen is distributed over a scale 
height about 10 times larger than the observed atomic hydrogen 
distribution, for the surface density to agree with the observed 
value, ii) it implies a scale length for the disk Rd <~ 10 kpc, 
about 3 times larger than the value deduced from star counts, 
for the circular rotation velocity to be constant up to ~ 30 kpc, 
iii) it does not explain the observed microlensing events toward 
the LMC. In case such clouds are distributed in the halo (Ger- 
hard & Silk, 1996), point i) and iii) remain valid. An exten- 
sion of this model, namely a highly flattened (q ~ 0.1) halo, 
although not completely ruled out, is at odd with various ob- 
servational constraints (see §2). 

Two models (and possibly a whole range of models in- 
between these extremes) fulfill all the afore-mentioned obser- 
vational contraints. The various components of these two mod- 
els are given in Table I. Both include: 

• a central, elongated bulge, with Mbui go ~ 2 x 10 10 M Q 
from kinematic considerations. The analysis of the first year of 
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MACHO observations towards the bulge (Alcock et al., 1997) 
seems to exclude a significant population of brown dwarfs in 
this region (Paper I) although better statistics are needed to re- 
ally nail down this issue, 

• a thin double-exponential disk with a scale length Ri ~ 3 
kpc and a scale height hi ~ 320 pc (Bahcall & Soneira, 1980) 
(or equivalently a sech 2 x exp model (Gould et al., 1997)), + 
a thick sech 2 (\z\/h2) disk with the same scale length, a scale 
height /12 ~ 640 kpc and a local normalization ~ 20% (Gould 
et al., 1997). The mass of the disk reads: 



M. 



disk 



2nZ Q R 2 d e R ° /Rd 



(16) 



where S Q is the local mass density, which is different in the 
two models. 

• a flattened (q ~ 0.6) DeVaucouleurs spheroid with a mass 
Msph ~ 1-3 x 10 10 M Q , i.e. a negligible contribution to the 
Galactic mass budget (Chabrier & Mera, 1997a), 

• a flattened dark halo with the density profile ( |i"2| ) 
yields a mass: 



which 



Mh 



1.36 x 10 



R, 

12 Rh 



Arctan^-) 

Rr 



23 kpc 



100 



M 



0! 



(17) 



23 kpc 

which is nearly independent of the value of the core radius R c . 

In both models, a significant contribution of low-mass stars 
and brown dwarfs to the dark halo mass is excluded from star 
count observations and microlensing experiments (Chabrier & 
Mera, 1997a; Mera, Chabrier & Schaeffer, 1997). Although we 
can argue that neither the EROS nor the MACHO experiments 
are optimized for short-duration (~ days) observations!, a S 'S~ 
nificant population of low-mass (brown dwarfs) objects in the 
halo would produce many short-time events: N = T x E x e, 
where E = 1.8 x 10 7 star.year is the exposure of the MA- 
CHO collaboration, e is their mean detection efficiency and 
r ~ 1.6 x 10~ 6 (m/M Q ) -1 / 2 event/year is the event rate 
for the standard halo model. A Dirac MF in the BD domain 
m ^ O.1M implies about 16 events with t e ^ 20 days, 
which is excluded at 99% confidence level, since only one has 
been observed (with t e w 19.4 days, Alcock et al., 1996b). 
Planetary-like objects (m < 10~ 2 Mq) are also excluded as a 
significant mass-contribution (<; 10%) at 95% confidence level 
by the EROS CCD observations (Renault et al., 1997) and the 
MACHO "spike" observations (Alcock et al., 1996a). In both 
models also, the dark halo contains a significant population of 
compact objects, which are responsible for the microlensing 
observations towards the LMC and amount to a local density 
p ~ 3 x 10~ 3 M Q pc~ 3 . As mentioned below, this condition 
is compulsory for the second (modified halo) model. Although 
uncomfortable, the halo WD solution represents the "least un- 
likely" scenario and cannot be excluded at the time these lines 
are written. 



4 Although the EROS CCD experiment was optimized for short 
timescales, the exposure, ~ 80000 stars x a few months, is small. 



Table 1. Main characteristic of the two consistent Galactic 
models outlined in the Conclusion. The dark halo is made of 
white dwarfs in both models, but the standard model has a dom- 
inant non baryonic component. 





Std halo 


heavy disk 


unit 


Disk: 








Scale length 


3.2 


3.2 


kpc 


Spfllp tlPlotlf 

iJ L lI 111 I L_ 111 


320 


320 


pv 


rtl 1 

h'ha 


< 015 


~ 0.03 


iviq 


Mass 


3 x 10 10 


5 x 10 10 


M^ 


Ttiiplc DisV' 

L 1 1 1 V. IV 1 ' 1 








Scale length 


3.2 


3.2 


kpc 


Spalp hpidtit 
tin. in. ii 


640 


640 


nc 
V- 


Mass 


6 x 10 9 


6 x 10 9 


Mn 


Bulge: 








Mass 


1.2 x 10 10 


1.2 x 10 10 


M Q 


Bar 


yes 


yes 




Brown dwarfs 


none 


none 




Spheroid: 








Density 


De Vaucouleurs 


De Vaucouleurs 




Oblateness 


q ~ 0.6 


q ~ 0.6 




Mass 


1.3 x 10 10 


1.3 x 10 10 


M 


Halo: 








Core radius 


5 


14 


kpc 


Composition 


< 50% WD 


up to 100% WD 




Mass 


~ 10 12 


~10 12 


Mq 




220 


220 


km.s -1 


Total column density: 








<E) at Rq 


~ 55 


~ 75 


Mq pc" 2 



The very difference between these two models resides in 
the location of the Galactic dark matter. The first model is 
simply what we call the "standard" model with a possible 
(~ 10±10 M© pc~ 2 ) amount of hidden mass under the form of 
brown dwarfs in the disk or in the bulge. In this case, nearly all 
dark matter is in the dark halo and is predominantly (> 50%) 
of non-baryonic nature. In the second case, the Galactic disk 
must necessarily include a significant brown dwarf popula- 
tion. The analysis conducted in Paper I (§4.1.2) yields a brown 



dwarf contribution £q bd ~ 20 — 3OM pc , comparable 
to the stellar contribution, and thus a local surface density 
S Q ~ 60 - 70M Q pc- 2 , since E 0vis = 43 ± 5M pc- 2 
(Eq. (15) of Paper I), within the la limits of the observed dy- 
namical value (Flynn & Fuchs, 1994). The average disk surface 
density at the solar radius is in this model ~ 5 — 20% larger 
than the dynamical value determined in the solar neighborhood. 
This may be due to the presence of density waves. Most of the 
dark matter is still in the dark halo but its contribution to the lo- 
cal dynamical density amounts only to~ 3 x 10~ 3 M Q pc~ 3 , 
~ 30% of the standard halo contribution. This (heavy disk) 
model necessarily needs the presence of white dwarfs (or other 
remnants) in the halo to be confirmed, these latter providing 
most of the missing mass. In this case, the missing mass in 
spiral galaxies is entirely baryonic and the dark baryons con- 
sist essentially of brown dwarfs in the disks and white dwarfs 
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in the halos. This model, however, must be considered as an 
extreme case. But it shows the possibility of a heavy Galactic 
disk, with a substantial fraction of brown dwarfs, which yields 
a smaller local normalization for the halo than assumed in the 
usual standard disk+halo model. This leads to a predominantly, 
rather than "entirely" baryonic model, a more likely solution. 
This model, however, requires a rather large surface density in 
the solar neighborhood and a relatively steep (a ~ 2) mass 
function in the disk near the hydrogen-burning limit, which ex- 
tends apparently (within the present uncertainties) to smaller 
masses than in the bulge (see Paper I, §4.1.1 and 4. 1 .2). 

Therefore, the key issue to determine the correct Galactic 
mass-model, based on the present analysis, is the determina- 
tion of i) the brown dwarf density in the disk/bulge and ii) the 
white dwarf (or other candidates) density in the dark halo. On- 
going microlensing experiments (EROS II, MACHO), infrared 
surveys and large field high-proper motion surveys might soon 
nail down this issue. Analysis of these observations in the con- 
text of the present models will yield immediatly the contri- 
butions of the various baryonic components to the different 
parts of the Galaxy and eventually the characterization of the 
(at least baryonic) dark matter in the Galaxy. Meanwhile, the 
DENIS and 2MASS projects should allow a precise determina- 
tion of the disk stellar mass function down to (perhaps below) 
the hydrogen-burning limit, allowing more robust extrapolation 
into the brown dwarf domain. Unfortunately, observations at 
the time these lines are written do not allow to disentangle be- 
tween these two models. Various uncertainties in present-day 
observations (scale heights, scale lengths, u rot , £©, ii©, mi- 
crolensing observations) remain within the error bars allowed 
by the present models. 

We have shown in the present paper that both types of mod- 
els yield the observed, constant circular rotation curve up to 
large galactocentric distances. In particular a truncated disk 
yields an exponential density profile. In the second model, the 
disk contributes to the rotation curve up to its edge, and the 
dark halo beyond this limit. Both the disk and the halo yield a 
flat rotation curve and there is no need anymore to invoke any 
peculiar conspiracy. We stress that, in virtue of the Occam's 
razor principle, we have considered in the present calculations 
the most natural hypothesis for the nature and the location of 
dark matter. We have assumed in particular that the dark mat- 
ter is distributed homogeneously, as visible matter. This yields 
about constant, but different, mass-to-light ratios in the disk 
and in the halo. We have not considered more exotic scenar- 
ios involving clustered dark matter in the halo (see e.g. Kerins, 
1997). We believe the present calculations to provide new in- 
sight on the possible nature and distribution of dark matter in 
galaxies. Comparisons of future observations with the present 
models should help answering this long-standing problem. 



Acknowledgments: We have benefited all along this work 
of numerous discussions with various people, more particularly 
with the MACHO and the EROS groups and with D. Pfen- 
niger and F. Combes. We also acknowledge useful conversa- 



tions with O. Gerhard and with the referee, A. Gould. It is a 
pleasure to thank these persons for their valuable help. 

Note added in proof. After the present calculations were 
completed, a new determination of the local dynamical den- 
sity has been derived from the Hipparcos data (Creze et al., 
fcstro-ph/9709022| ). This new determination leaves essentially 
no room for any dark component in the disk and thus, if con- 
firmed, seems to exclude the modified-halo model (§5) and to 
confirm the standard-halo model (§3). 
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